Solitons in a composite multiferroic chain 
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In a coupled ferroelectric/ferromagnetic system, i.e. a composite multiferroic, the propagation of 
magnetic or ferroelectric excitations across the whole structure is a key issue for applications. Of a 
special interest is the magnetic or ferroelectric soliton dynamics across the ferroelectric-ferromagnetic 
interface, particularly when the soliton carrier frequency is in the band of the ferroelectric and in 
the band gap of the ferromagnet. For a proper choice of system parameters, we find that there is 
a threshold amplitude above which the interface becomes transparent and a band gap ferroelectric 
soliton penetrates the ferromagnetic array. Below that threshold the soliton is fully reflected. Slightly 
below the transmission threshold, the addition of noise may lead to soliton transmission, provided 
that the noise level is not too low or too high, as found in stochastic resonance. These findings 
represent an important step towards the application of ferroelectric and/or ferromagnetic soliton- 
based logic. 

PACS numbers: 85.80.Jm, 75.78.-n, 77.80.Fm 



Introduction - The revival of research on multiferroics 
(MP), i.e. systems with coupled ferroic (magnetic, elec- 
tric, or elastic) orders [l|-[3| is fueled by recent advances in 
synthesizing composite ferroelectric (FE)/ferromagnetic 
(FM) nano and multi layer structures that show a sub- 
stantially larger multiferroic coupling strength [l|-[5| as 
compared to bulk matter, i.e. single-phase multiferroics 
[gI, Q such as Cr203[8[. MF are important for addressing 
fundamental questions regarding the cross correlation be- 
tween electronic correlation, symmetry, magnetism, and 
polarization but also hold the promise for qualitatively 
new device concepts that exploit magnetoelectric (ME) 
coupling to steer magnetism (ferroelectricity) via electric 
(magnetic) fields. Potential applications are wide ranging 
from sensorics and magnetoelectric spintronics to envi- 
ronmentally friendly devices with ultra low heat dissipa- 
tion (ol-llH . Thereby, a key issue is how efficient magnetic 
or ferroelectric information, i.e. an initial excitation, is 
transmitted in a system with a MF coupling. For in- 
stance, in a two-phase or composite MF [6f,0,ll3 such as 
BaTi03/CoFe204 
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^ „, PbZri-xTixOs/ferrites PJ, |15|, 
BaTiOa/Fe [16|, PbTiOs/Fe EES or BaTiOs/Ni the 
MF coupling is strongest at the FE/FM interface whereas 
away from it the EE or FM order is only marginally af- 
fected. We expect, thus that a ferroelectric signal trig- 
gered by an electric field in the FE part may or may not 
be converted into a magnetic signal depending on the dy- 
namics taking place at the interface. How this transport 
of information depends on the properties of the system 
is rarely studied and will be addressed in this work. The 
outcome of such a study would not only uncover the con- 
ditions for optimal signal handling but also holds the 
potential for new insights into the mutliferroic coupling 
retrieved by tracing the signal dynamics. We will focus 



on nonlinear waves forming solitons that are known, in 
the pure FM or FE phase, to have a series of applica- 
tions in magnetic logic, microwave signal processing, and 
spin electronic devices. Advantageous in this respect is 
that solitons (i.e. localized state of a large number of 
elementary excitations, e.g. magnons) are very robust 
and have a particle-like nature [19]. Hence, logic oper- 
ation based on magnetic soliton attracted detailed con- 
siderations [20I, HH. Multiferroics offer new fascinating 
mechanisms for the soliton dynamics. For example, due 
to discreteness and/or nonlinearity of the system, it may 
happen that the frequency of the large amplitude exci- 
tation may fall within the gap of the linear oscillations 
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FIG. 1. (color online) Schematics of a chain consisting of fer- 
roelectric and ferromagnetic parts coupled at the interface. 
Lower panel shows a particular choice of frequency for which 
in the ferroelectric a conventional envelope soliton is formed, 
while in ferromagnetic part bandgap soliton develops with a 
nonlinear frequency shift proportional to (B is an ampli- 
tude of magnetic gap-soliton). The chosen mutual alignment 
of the ferroelectric polarization and the magnetization at the 
interface resembles the realized ferroelectric (BaTiOs) tunnel 
junction with ferromagnetic (Fe) electrodes [22l] . 
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spectrum, as illustrated in FiglH Therefore, the energy of 
the excitation would not spread over the lattice. As well 
established from detailed studies on intrinsic localized 
modes, e.g. (23|-[3Q|, we know that in spite of the local- 
ized energy profile, such modes may move along the whole 
chain. This means that excitations created in the ferro- 
electric part via an electric field can be transmitted to the 
magnetic part and move there. Moreover, as it will be 
shown below the creation of these band-gap localizations 
could be enhanced by the noise manifesting thus some 
analogy with stochastic resonance phenomenon |3lL l32|. 

Model- For our purposes a large ME coupling is nec- 
essary. In this regard, new fabrication methods [15| for 
so-called two-phase or composite multiferroics [Sl, lD, 12 | 
as well as the realization of ferroelectric wires [33 1 
are encouraging. Examples of composite multiferroics 
BaTi03/CoFe204 [i3 or PbZri.xTixOa/ferrites [TBI are 
still popular. Major research is focused on BaTiOs/Fe 
[Tel, PbTiOs/Fe Q^,^ or BaTiOs/Ni composite multi- 
ferroics, to name but a few, since their bulk parameters 
are very well known (Ref. [3J for BaTiOs and Ref. [35| 
for Fe or Ni) as well as the misfit of lattices is relatively 
low Relatively high ME constants 17, 18| were pre- 
dicted for these materials at room temperatures. A possi- 
ble mechanism for ME coupling at the FE/FM interface, 
that is also assumed here, is based on screening effects 
(36j : In the polarized FE layer surface charges are cre- 
ated which are compensated by screening charges in the 
metallic FM electrodes when attached to the FE layer. 
In this heterostructure the accumulated screening charge 
density on the FM side of the interface is spin-polarized 
and couples then to the magnetization in FM. Techni- 
cally, we describe the bulk unstrained BaTiOs by the 
Ginzburg-Landau-Devonshire potential FGLD [34] . For 
the discretized FE polarization (P^) in a coarse-grained 
approach, the form of the FGLD potential for a gen- 
eral phase and arbitrary temperatures is quite involved 
(37| . However, at room temperature BaTiOa-crystal has 
an axis along which the polarization switches (tetrago- 
nal phase). Consequently, the form of the FGLD poten- 
tial reduces to the one dimensional biquadratic poten- 
tial. For the description of the magnetization dynamics 
in iron (Fe) {Sk, discretized and normalized to the satu- 
ration value of the coarse-grained magnetization vector), 
we employ the classical Heisenberg model. With the aim 
to explore the feasibility of conversion of the electric ex- 
citation formed in FE part into a localized spin magnetic 
excitation in FM part we employ thus the multiferroic 
model (cf. Fig. [1]) 



where Hp is the Hamiltonian of FE part for the mul- 
tiferroic system of A^-interacting FE dipole [H, [s^ Pn- 
Note, that at room temperature, the polarization vec- 
tor has only one nonzero component and we choose 
the polarization vector to be directed along the x axis 
Pn = (Pn, 0,0), n = 1,...,A'. ao is a kinetic coeffi- 
cient, ai^2 define the FGLD potential constants and is 
the nearest neighbor coupling constant. Hs is the classi- 
cal Hamiltonian of the ferromagnetic part of the system 
[ssj . where J is the nearest neighbor exchange coupHng 
in the FM part. The coarse-grained FM order param- 
eters are 5^, k = A' + 1,...,M and D is the uniaxial 
anisotropy constant. Interface effects between the spin 
and the FE dipole systems are described by the dipole- 
spin interaction Hamiltonian Vsp- Since we assume the 
ME coupling to be due to screening and the screening 
constant in metals is on the scale of the lattice constant 
we couple only the (coarse grained) dipole at the interface 
(indexed N) to the interfacial spin (with index 1). 

In our numerical simulations we will operate with di- 
mensionless quantities upon introducing Pn = Pn/Pmaxi 
Sk = Sk/S and defining a dimensionless time as t ^ ujot 
{ujo = y^tv/ao ~ 10"^^ Hz). The equations governing the 
time evolution of the dipoles and the spins (except for 
the sites near the interface) read 



dt2 



-apn - f3pl + {pn-l - 2pn + Pn+l) (2) 



dt 



±iJ [4 (4-1 + 4+i) - s| (4-1 + 4+i)] ± 



(3) 



where n ^ N and k ^ A'+l. We introduced the following 
dimensionless constants a = ai//^, /3 = a2P^axl^^ J — 
JiS/uoq and D = J2S/UOQ. For the dipole pn and the spin 
si at the interface the relation applies 



d PN 



-apN 



(3p^ 



N 



^{PN-i-2pN^gssf) (4) 



ds 
dt 



^ = ±iJ [sfs^ - sfsf] ± i [2Dsfsl - QpPNsl] . 



Here = s% ± isl, Qs = gS/{nPmax) and = 
gPmax/ {Sujq). The evolution according to Eqs. (j2]|4]) 
proceeds under the constraint (s^)^ + {^\)^ + {^kf' ~ ^• 
For the derivation of weakly nonlinear envelope solu- 
tions from eqs. (j2j) and (j3j) one can utilize the reductive 
perturbation theory developed in Ref. [H, |4o| and ob- 
tain the solutions for the dipoles and spins separately in 
the following form (a detailed derivation is provided as 
supplementary materials to this paper): 



H = Hp + Hs + V^sp, 



(1) 



Hs= i-JiSkSk+i - J2 {Sif). Vsv = -gPNSf, 

k=N-\-l ^ ^ 



Pn 



A cos [ujpt — QpU -\- Sujpt] 

cosh [(n - Vpt)/Ap] 



~ cosh [(A: - Vst)/As] 
(5) 

where A and B are the amplitudes of the dipolar and the 
magnetic envelope solitons, respectively; Up and ujs are 
the frequencies of the linear excitations which obey the 
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following dispersion relations 



cjp = y a + 2 (1 - cos^p), cjg = 2 [I) + J(l - cos^g)] , 

(6) 

qp and are the carrier wave numbers of the dipolar and 
the spin excitations; Vp = ^mqp/ujp and Vs = 2 J sin 
are the group velocities of the corresponding solitons. 
The width of the dipolar and spin solitons are 

The nonlinear frequency shifts are defined as 

Suj,=A'r^, Sco.s = -B''^. (8) 

At first we note that for the wave packet transmission 
the matching condition between the frequencies has to 
be fulfilled 

ijjp + 6ujp = + Sujs' (9) 

In addition, for an efficient transmission of the soliton 
from FE into FM part the widths of the solitons should 
be the same in both parts, i.e. Ap = with the restric- 
tion B < QpA. If one excites the soliton with a carrier 
frequency uj which is located within the band of both 
the dipolar and the spin wave spectrum, then the soli- 
ton will safely penetrate from the FE to the FM part, 
but some portion of the energy will be reffected from 
the interface. Changing the soliton amplitude one can 
manipulate the ratio between the transmitted and re- 
ffected parts. In addition, the transmission is very sen- 
sitive to the coupling constant g between FE and FM 
parts. We have investigated this dependence by ffxing 
the values for the dimensionless parameters as follows: 
a = 0.2; f3 = 0.1; J = 1; D = 0.6 and varying the 
coupling constant gg. We assume for simplicity gp = gs 
(in general these constants differ depending on the sam- 
ples material but present no obstacles for theory). 

Numerical Results.- Graphs a), b), c) of Fig. [2] show 
the localized energy evolution along the lattice for differ- 
ent values of coupling constant gp. In graph d) the de- 
pendence of transmitted energy on the coupling constant 
is displayed evidencing that the transmission is maxi- 
mal when gp lays in between the spin coupling constant 
J = 0.5 and dipolar coupling constant (in reduced units 
it is equal to 1). An insight into the spins alignment as 
well as the topology of the excitation and its propagation 
in the chain's FM part at different time moments is dis- 
played graphically in the supplementary material to this 
paper. 

Further interesting effects arise when a band soliton 
forms having a carrier frequency uj located in the band 
of the dipolar spectrum (see bottom panel of Fig. [1]) and 
slightly below the zone boundary uJs{qs = 0) of spin wave 
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FIG. 2. (color online) Insets a), b) and c) show the depen- 
dence of the local energy (for dipoles it is given by the local 
values of Hp and for spins it is obtained from the local values 
of the expression Hs + S^[Ji + J2]) on the site number and 
the simulation times as dictated by Eqs. (|2])-(|3]) for the set 
of parameters given in the text and for the different values of 
coupling constant gp that are indicated in the graphs {gp oc g 
and g is the interfance ME coupling strength). The dipolar 
soliton carrier wave number is chosen as Qp = OAtt and the 
dipolar soliton amplitude is chosen as A = 0.2. Dipoles and 
spins are occupying the sites n — 1 ... 150 and k = 151 .. . 300, 
respectively, d) The dependence of the relative energy trans- 
ferred to FM part (i.e. the ratio of the energy in FM part on 
the whole injected energy) on the coupling constant strength 
9p- 

spectrum. Then for small amplitudes the dipolar soliton 
is totally reffected from the interface because it does not 
resonate with any mode in the spin array. However with 
increasing the amplitude, there is a threshold value Acr 
(due to the nonlinear frequency shift) above which the 
soliton starts to transmit towards the FM part of the 
multiferroic chain forming thus a magnetic gap soliton. 
Using Eqs. (|8|) and (|9|) and assuming B = A one can 
infer the relation deffning this threshold amplitude to be 

..i,.=o)-.=(^^+^yi. (10) 

Based on this observation, we proceed with the simu- 
lations according to Eqs. (|2])-(j4|) with the set of param- 
eters given above and choosing gp = g^ = I and starting 
at t = with a soliton in the form of the ffrst expression 
in Eq. (j5j) with a carrier wave number qp = 0.37. For 
such a wave number the corresponding linear frequency 
is a; = 1.1856 thus it is located in the band gap of the 
spin wave spectrum and no soliton transmission happens 
in the case of small amplitudes as it is seen from graph a) 
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FIG. 3. (color online) The same as in graphs a)-c) of Fig. [2]for 
the set of parameters given in the text and Qp = Qs = I. The 
dipolar soliton carrier wave number is chosen as qp — O.STtt 
and the dipolar soliton amplitude values are indicated in the 
graphs. Color scale as in Fig. [21 

of Fig. m According to the relation (p!Q|) we can calculate 
the threshold amplitude for which a soliton transmission 
emerges and find Acr = 0.22. In the numerical results 
the transmission occurs for the incident soliton ampli- 
tudes A > 0.27. This discrepancy is explained by the fact 
that the solitons amplitudes in different parts of multi- 
ferroic structure does not exactly coincide. In the graph 
b) of Fig. [2] we display the dynamics for a larger soli- 
ton amplitude, i.e. A = 0.33 and find that an envelope 
soliton is formed in the ferromagnetic part as well. With 
further increasing the soliton amplitude, the transmitted 
soliton takes over almost all the energy of the incident 
one (graph c) of the same figure). 

keeping the amplitude of the incident soliton slightly 
below the threshold (in our simulations we choose A = 
0.265) it is obvious that even a small perturbation may 
cause transmission to FM part. Thus we add a term 
f{t)Si to the Hamiltonian (j2j) describing the action of 
a random magnetic field at the interface spins. f{t) is 
uncorrelated in time and randomly distributed in the in- 
terval [—/,/]. Then for small random fields / = 0.05 
the picture is almost the same as in the case of a zero 
noise (please compare upper graphs of Figs. [3]and|l]). 




FIG. 4. (color online) The same as in Fig. [3] (including the 
color scale), however with a fixed soliton amplitude A = 0.265 
and different noise levels indicated in the graphs. 



Increasing the noise strength to a moderate level, energy 
transmission in FM part takes place (see graph b) of Fig. 
HI). This stochastic resonance like behavior is displayed 
in graph c) of the Fig. HJ 

Summary- As evidenced by analytical and numerical 
results, in a two-phase multiferroic the magnetoelectric 
coupling at the interface renders possible the conversion 
of an initial ferrolectric soliton into a ferromagnetic sig- 
nal, opening thus the door for FE and/or FM soliton- 
based logic in multiferroics. As an essential step in this 
direction we identified the conditions under which a FE 
signal is converted into a FM one and clarified the role 
of noise. 
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SUPPLEMENTARY MATERIAL 



As an additional detail to figure 2 of the main text, Fig l5] provides an insight into the real space spin alignments 
and their time evolution for the same parameters discussed for figure 2 of the main text. 



t=150 




FIG. 5. (color online) The real space spin configuration in the FM part of the chain at certain time moments (t — 150 and 
t = 300 in normalized units) displayed for the for the case depicted in Fig. 2(c) in the main text, sf^ + sfy at the sites i of 
the FM chain is also shown. 
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MATHEMATICAL DETAILS 



Furthermore we make the reduction of Pn+i and Pn-i 



A solution method for Eq.(2) of the main text is 
needed. This equation reads 



-apn - + {Pn-l - '^Pn + Pn+l) (H) 



which describes the FE dynamics that we capture via 
the reductive perturbation method [39| that builds on a 
multiple scale expansion in the following form 



oo oo 



,im{ojt—qn) 



(12) 



a=l m— — oo 



where ^i^^ = ( ^^2, ) and the slow variables are 



^ = e{n-Vt), T = e^t. 



(13) 



We insert (p!2]) into the initial equation ([TT]) and keep 
the expansion terms up to the third order over e. For 
this the following derivative over time is needed 

dt dC"" dr ' ^ ^ 



oo oo 



Pn±i = ^e'^ ^ M/(;^)(C±£,r)e^™("*-«"T«), (15) 



a—1 m— — oo 



and perform the expansion of ((f ± £,r) in terms of 
small e. We recall also that only terms up to third order 
over e need to be kept. Then we obtain the following 
expressions for Pn and Pn±i 



m— — oo ^ ^ 



^im{ojt — qn) 



Pn±l 



m— — oo 



im{ojt — qn^q) 



(16) 



Inserting these into Eq. (fTT]) and taking into account the 
relation (fT4|) for the time derivative we obtain to a first 
approximation with respect to e the following relation 
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E 



^im(ujt-qn) r 2^2 



[w^m^ - a + (e™« + e-™« - 2)]*(^) = 0. 



(17) 



Let us assume that is nonzero. From this the dispersion relation fohows straightforwardly 



u;^ = a + 2(1 — cosg) 



(18) 



and all other with m 7^ ±1 are exactly zero. 
In the second approximation one has the relation 

00 00 



-tmq %mq 



)] 



(1) 



m= — cx) 



- = 0. 

(19) 

From the dispersion ([18]) and the fact that all with m 7^ ±1 are zero it follows that with m 7^ ±1 are zero 
as well and one has the second relation defining velocity V 

V = sing/cj 

Finally in the third approximation with regard to e we find 



(20) 



^t-qn) I [^2^2 _ ^ _^ (-gimg _^ g-img _ 2)] ^(^) + 1 (_2y2 _^ g 

(2iwmy + e-™« - e'""*) - 2imw^|^| = 



imq _|_ p — imq 



(21) 



where the nonlinear term in the third approximation reads 



{pp. 



n) 



3^ I ^^(1)^^ g3i(u;t-gn) ^ ^ ^ g-3i(a;t-gn) ^3 ^(1) ^ ^{^) ^i{ujt-qn) ^ ^*(^)e-*(^^-?^)^ j(22) 



Collecting the terms proportional to e*^^^ (meaning that we choose the harmonics with m = 1) and taking into 
account the relations (p!8|) and (|2Q|) we obtain 



(1) 



One soliton solution of the nonlinear Schrodinger equation ([23]) can be written in the form 



(23) 



= exp 



'.3(3 ■ 
4uj 



cosh 



(24) 



and thus in the first approximation over e we find the following expression for the coupled dipoles according to the 
expansion ([T2]) 



Pn = 2e cos 



ut — pn-\- e — t 
4uj 



cosh 



e{n- Vt) 



b:^ - o? - 4a 



(25) 



or assuming, that the amplitude is defined as A = 2^ 



Pn = A COS 



,2 3/? 



ujt — pn-\- A 1 



cosh 



A{n-Vt)^ 



2 (a;4 - a2 - 4a) 



(26) 



which is the same as the first expression in Eq.(5) with the definitions Eqs.(7) and (8) of the main text. 
For the spins we rewrite the equation Eq.(3) of the main text 



ds 
~dt 



^ = ±zJ [4 {sl_, + - si {sti + 4+i)] ± ^^Ds^s 



k 



(27) 



where the z spin component is defined from the relation 5^5^ + ("^D^ = 1- We proceed again with the similar 
expansion as in the dipolar case 



4 = E^" E *L"n?,T)e^ 



im{ojt — qk) 



1 



Sjl s 



k '^fc 

2 ' 



T eh 



(28) 



a=l m— — oo 

Substituting these into the initial equation (|27|) we find to a first order in e the following relation 



e e 

m— — oo 



im{ojt — qk) 



[muj - 2L> + 2 J (cos mq - 1)] = 0, 



assuming that ^^^^ is nonzero, one obtains the dispersion relation 

a; = 2L> + 2J(l-cos(7) 



(29) 



(30) 



and the condition that all with m ^ 1 are equal to zero. In the second order we choose the harmonics e'^^^^ 
(i.e. the harmonics with m = 1) and find 



2i:>-2J(l -cosg)]^f^ +i(V- Jsing) 



(1) 



0. 



Taking into account the dispersion relation (|3Q]) one arrives at the relation for the soliton velocity 

V = 2 J sing. 

In the third order approximation, we collect the harmonics with m = 1 and obtain 



(31) 



(32) 



[uj-2D- 2J(1 - cosg)] -i{V - 2 J sing) 



J cosq 



Thus, taking into account (pO|) and ()32|) the nonhnear equation has the form 



[D + J(l-cos</)] 



(1) 



J cosq 



dr ^ ^^^^ 2 



(1) 



0. 



In full analogy with the dipolar case, one soliton solution of this equation has the form 

e{k - Vt) 



^^^^ exp 



cosh 



4 J COS g 



= e exp 



cosh 



4: J cosq 



(33) 



(34) 



Defining the soliton amplitude B = e we recover the second expression appearing in the main text under Eq.(5) with 
definitions Eqs.(7,8). 

I 



